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Abstract 

Casorati determinant solution to the non-autonomous discrete KdV equation is constructed by using the bilinear 
formalism. We present three different bilinear formulations which have different origins. 



1 Introduction 

In this article, we consider the following partial difference equation 

1 + ^],"' + (1) 

am bn+\l \flm+i b„ j " \a„ b„j v"^ \fl„,+i b„+il v"^l^l 

where m, n are the discrete independent variables, V" is the dependent variable on the lattice site (m, n), and fl„„ b„ 
are arbitrary functions of m and n, respectively. Recently eq.© has been derived by Matsuura|l| as the equation of 
motion of discrete curves on the centro-affine plane. In this context, v"' is related to the curvature and it is essential 
that a,„ and b„ depend on m and «, respectively. For analyzing the motion of discrete curve, constructing solutions of 
eq.([T]) explicitly is an interesting subject. We call eq.([T]i the non-autonomous discrete KdV equation, for if a,„ and b„ 
are constants, e.g. a,„ = a, b„ - b, eq.© reduces to the discrete KdV equation f2l[3l 



or 



m+l ^ 

a + b 



ni+l 



(3) 



Since the constants a and b correspond to the lattice intervals of m and «, respectively, eq.([T]i can be also regarded as 
the discrete KdV equation on inhomogeneous lattice. 

The non-autonomous version of discrete integrable systems on two-dimensional lattice have not been investigated 
well, although those on three-dimensional lattice, such as the Hirota-Miwa (discrete KP) equation or the discrete two- 
dimensional Toda lattice equation, have been studied well together with their solutions|4, 5 , 6, 7, 8|. Let us explain 
a reason taking eq.([T]i as an example. The autonomous version eq.© can be transformed into so-called the bilinear 
equation by suitable dependent variable transformation. The bilinear equation is regarded as a reduction of the Hirota- 
Miwa equation|9, 10|, which is well-known to admit various types of exact solutions, such as soliton solutions|3, 9|, 
rational solutions expressible in terms of the Schur functions lflOl ITTII . or periodic solutions that are written in terms 
of the Riemann theta functions! 12] . Therefore one can obtain solutions to the discrete KdV equation (|2| by applying 
the reduction procedure to those for the Hirota-Miwa equation. Now, the Hirota-Miwa equation and its solutions 
can be generalized to non-autonomous case in a straightforward manner. However, it is shown that one cannot apply 
the reduction procedure to the non-autonomous Hirota-Miwa equation consistently. Moreover, eq. ^ cannot be 
put into bilinear equation by the procedure similar to the autonomous case because of the non-autonomous property. 
Therefore it was not clear how to construct solutions to the two-dimensional non-autonomous discrete integrable 
systems systematically. 

In this article, we construct the Casorati determinant solution to the non-autonomous discrete KdV equation ([T]) by 
using the bilinear formalism. We present three different bilinearizations: The first one can be derived by the reduction 
of non-autonomous discrete KP hierarchy with a new technique. The second one is the bilinearization obtained by 
introducing certain auxiliary r function which has a similar structure to the ones that appeared in the study of Rj and 
Rii biorthogonal rational functions lfTSl [141 . The third bilinearization is through the use of non-autonomous potential 
discrete KdV equation. 



This article is organized as follows. In Section 2 we review the bilinearization of the discrete KdV equation 
(|2]l, and discuss briefly why the similar calculation fails for the non-autonomous case. In Section 3 we discuss the 
bilinearizations of eq.([TJ and construct the Casorati determinant solution. Finally, concluding remarks are given in 
Section 4. 



2 Bilinearization of the discrete KdV equation 

The discrete KdV equation ^ can be transformed to the bilinear equation 
by the dependent variable transformation 

_m 1 

m = "+1 " (5) 

In fact, substituting eq.® into eq.© we have 

4- Mt™ -m+\m+2 ( 1 , M m+2 m m+l 
\ - ^jClC'^r' - - (6) 

Interchanging the second term of the left hand side and the first term of the right hand side, and dividing the both sides 

by C+V';r;f+' we get 

( i + 1 ) t'" t^V - ( i - f ) t'" , t"'V (- + \) Tt,V'+2 - ( i - i) r^+iT™:^ 

\a b ) «+l «-l \a b ) n-\ n+1 _ \ b ) n+1 «-l \ a o/ n~\ n+1 

_m_m+l _m+l _m+2 

^ n hi ^ n n 

Equation can be decoupled as 



since the right hand side of eq.O is obtained from the left hand side by shifting m to m + 1 . Here a(n) is an arbitrary 
function in n, which can be absorbed by suitable gauge transformation on r™. We obtain eq.© by choosing a(n) = | 
so that = 1 is a solution. 

The bilinear equation (|4]i can be obtained by applying the reduction to the Hirota-Miwa equation 

fli(fl2 - 03) t(Ii + 1, 12, It,)t{Ii, I2 + 1, Z3 + 1) + 02(03 - fli) t(Ii, I2 + 1, h)T{l\ + 1, 12, h + 1) 
+03(01 - 02) t(/i, ^2, h + IM^i + I,l2 + 1, ^3) = 0, (9) 

where oi, 02, 03 are arbitrary constants. In fact, imposing the condition 

t(/i + 1,Z2 + 1,/3) = t(Zi,/2,/3), (10) 

where ~ means the equivalence up to gauge transformation, using eq.dTOll to suppress the li dependence and putting 
ai = -02, eq.® yields 

- (a2 - Oi) Tik - l,l3)T(l2 + l,h + 1) + (02 + 03) t(/2 + 1,/3)t(Z2- 1,/3 + 1) 

-2o3r(Z2,/3 + l)T(/2,W = 0, 

which is equivalent to eq. (|4|i with I2 - n, Ij, - m, 02 = b, aj, - a and t(/2, /3) = tJJ'. 

Now let us consider the non-autonomous case. We show that neither direct bilinearization nor reduction from the 
non-autonomous Hirota-Miwa equation work successfully for this case. First, substituting eq.© into eq.([T]l and doing 



the same calculation as above, we arrive at the following equation 



+ 1 

n+l 



—m+l—m+2 ' 



(12) 



which cannot be decoupled into the bilinear equation because of n dependence of the coefficients. Therefore naive 
bilinearization fails for the non-autonomous case. 

Secondly, let us consider the reduction from the non-autonomous Hirota-Miwa equation||7l[8l 



ai(h)(a2(l2) - 03(^3)) T(li + 1, h, h)T{h,h + 1, /3 + 1) 
+a2(fc)(a3('3) -fli(^i)) Tduh + 1,/3)t(/i -i- 1,/2,/3 + 1) 

-Hfl3(Z3)(fli(Zi) -02(^2)) HluhJi + IM/l + l,/2 + 1,^3) = 0, 



(13) 



where fl,(/,) (/ - 1,2,3) are arbitrary functions. Imposing the condition (fTOl i on eq. JTSb and suppressing the h- 
dependence, we obtain two different bilinear equations 



fll(^l)(fl2(/2) - 03(^3)) Tih - 1, /3)t(/2 + 1, /3 + 1) 

+a2(l2)(a3(h) - aiih)) rih + 1, /3)t(/2 - 1, ^3 + 1) 

-l-fl3(/3)(fll(/l) - 02(^2)) Tik, h + l)T(l2, h) = 0, 

adh - l)(fl2(^2 - 1) - 03(^3)) Tih - l,h)T(l2 +l,h + 1) 
+02(^2 - l)(fl3(/3) - aiih - 1)) Tih + l,h)Tih - 1, ^3 + 1) 

+fl3(/3)(fll(/l - 1) - fl2(/2 - 1)) Tih, h + IMh, h) = 0. 

Since those two equations should be equivalent, the coefficients must satisfy 

aiihKaiih) - a^ih)) a2ih)iaiih) - a^ih)) 



(14) 



(15) 



a\ih - \)ia2ih - 1) - fl3(/3)) 02(^2 - l)(fl3(^3) - a\ih - 1)) 

a3(/3)(ai(^i) - ^2(^2)) 
fl3(^3)(fli('i - 1) - 02(^2 - 1))' 



(16) 



which yields 



03(^3) 



1 



1 



fli(Zl)fl2('2 - 1) fll(/l - 1)02(^2)^ 
1 1 \ / 1 1 



fli(Zi - 1) aiih)/ \fl2(^2 - 1) a2ih) 



(17) 



Since this should hold for any aih) for all /3, we deduce that aih) and aih) must be constants. This implies that it is 
not possible to impose the condition (fTOl i on the non-autonomous Hirota-Miwa equation ( fT3] ) consistently, unless it is 
reduced to the autonomous case. 



3 Bilinearizations of the non-autonomous discrete KdV equation 
3.1 Reduction from the discrete KP hierarchy 

The non-autonomous discrete KP hierarchy in the bilinear form is expressed as|l2l[3] 

1 fl,-,(/„) «,,(/,, )2 •■• «,,(/„)"'-' a,,ihr-^- 



1 ai^li^) ai^U^f 



1 aijhj aiJh„f 



ai.ikr-^ 



aiSkT-^ 



TiJi, 



ai„,iliX'' aiJkr-'Ti^T. 



= 0, 



(18) 



where {; 



,im] C {1, 



,n|, T4 andr,; = 1, 



, m) are given by 



-T T 

;i /^ 



■ T T 



■T,T, 



(19) 
(20) 



respectively, ay{ly) are arbitrary functions in ly for each v, and n and m are arbitrary integers satisfying n > m > 3. 
Here T, is the shift operator of Z, defined by 



r,T(/i, /2, ■ ■ ■ , /„) = THuk, ■■■ , U-\,k + 1, U+\, ■■■ , In)- 

The simplest equation in the hierarchy (m - 3) is the non-autonomous Hirota-Miwa equation 

ai(li)(aj(lj) - ak(lk)) t(// + 1, Ij, hMk, lj + l,k+l) 
+a j(lj)(ak(h) - Oidi)) T(li, Ij+l, k)T(k +\,lj,k + l) 
+ak(h){ai{li) - ajilj)) rih, Ij, k + ^Mk +\,lj+\,k)^Q, 



(21) 



(22) 



where {/, j,k] c {!,...,«} and we suppressed other independent variables. The Casorati determinant solution to the 
hierarchy can be written as 



t(Ix,--- ,l„) = 



<f^^\h,--- ,ln) <p'^^"(lu--- Jn) 



■■Jn) 'P\'^'\lu 
■■Jn) v'C'\h, 



Jn) 



(1 



where (p^f^ {r - I, 



. , N) satisfy the linear equations 



, ly, ■ ■ ■ , Ifl) 



ay(ly) 



X) 

Jn) 
Jn) 



Jn), 



(23) 



for V = 1 , 



, n. For example, the A'-soliton solution is obtamed by choosmg ^J. as 

n ly-\ n ly-l 



''^;\h, ■■■Jn) = arpl Y[ Wd + ay{i)p,) +f5,.q'r n + 



(24) 



(25) 



V=l i-iy 



V=l i=iy 



where ar,l3r, p,-, qi (f - 1, ■ • ■ , A^) are arbitrary constants. 

Let us consider the reduction to the non-autonomous discrete KdV equation. The key idea is to specialize some 
of the independent variables to be autonomous by choosing the lattice intervals as constants and use them as auxiliary 
variables. Then one can use the autonomous variables for the reduction procedure to get the non-autonomous discrete 
KdV equation. We consider the four independent variables k - liJ - I2, m = 1^, n - I4 with the lattice intervals being 
6 = ai(k), £ = 02(1), ai„ = a^im), b„ - 04(11), respectively. We note that 6 and e are constants, namely, k and I are 
autonomous variables. Then we have the following bilinear equations from eq.i 



6(a,„ - b„) T(k + lj,m, n)T(k, l,m+ l,n + 1) 
+a,„{b„ — S) T(k, l,m + l,n)T(k + \J,m,n + 1) 
+bn{6 — am) T(k, I, fn, n + l)T(k + lj,t7i + l,n) — 0, 

e{a„, — b„) T(k, I + 1, m, n)T(k, l,m+ 1, « + 1) 
+a„(b„ — e) T(k, l,m + \ ,n)T{k, I + l,m,n + 1) 
+b„{e - a„,) rik, I, m, n + \)T{k, I + \,m + \,n) - Q. 



We impose the condition 



T{k + \J+ l,m,n) ~ T(k, I, m, n). 
This is achieved by imposing the condition on ipf'' (r - 1 , . . . , A^) as 

(fi'-^^ik + IJ + l,m,n) ~ (fil^^ik, I, m, n). 



(26) 

(27) 

(28) 
(29) 



For the case of soliton solutions, ipf\k, I, m, n) is expressed as 



/, m, n) ^a,pl{\ + dp^fH + ep,-)' + a^p,) + bjp,) 

i=mo j=riQ 
m-l n-l 

+/3,-qlil + dqr)\l + eqri ]~[(1 + fl,-^,) ]~[(1 + hflr). 



In order to satisfy eq.(|29]l, one may take 
so that 



qr = -pr, 5 = -e, 



^}~Sk + 1, Z + 1, m, n) = (1 - e^p^^ ip^;\k, I, m, n), 

N 

T(k l,m, n) = - e^pl) T{k, I, m, n). 



Then, suppressing the A;-dependence by using eq.(|28l). the biUnear equations ( |26] l and ( |27] l are reduced to 

— e(am - bn) t(1, m, n)T(l + 1, m + 1, n + 1) + a^ibn + e) t(Z + 1, m + 1, n)r(Z, ot, n + 1) 

— Z7„(e + fl,„) t(/ + 1, m, « + 1)t(/, m + 1, «) = 0, 

e(fl,„ — bn) t(1 + 1, m, m)t(/, m + 1, n + 1) + Omib,, — e) r{l, m + I, n)T(l + \,m,n+ 1) 
+ b„{e — a,„) t(1, m,n + 1)t(Z + 1, m + 1, n) = 0, 

respectively. By putting 

tJ" = t(/, m, «), o-"' = t(/ + 1, m, n), 
the above bihnear equations are rewritten as 

- e(fl,„ - bn) + ajb„ + e) Ci'^r' - + ««) ^r'<+i = 0, 



0, 



(30) 



(31) 



(32) 

(33) 
(34) 



respectively. Equations ( |33] ) and ( |34| ) can be regarded as a bilinearization of the non-autonomous discrete KdV equa- 
tion ([T]i. In fact, introducing the variables ^F"' and v™ by 



—tn _m+l 

" _/» 1 ' 
'n+1 



we obtain 



1 1\ 1 



1 1 



1 1 



r" + 
b n I " 



e bn] U„, e 



(35a) 
(35b) 

(36) 
(37) 



which are regarded as the auxiliary linear problem for the non-autonomous discrete KdV equation. Eliminating 
by considering the compatibility condition we obtain the non-autonomous discrete KdV equation ([TJ. The A^-soUton 
solution is given by 



ip\ (m, n) ifi^ (in, n) 
^2 ym, n) (in, n) 



(.S)/ ^ (s+1)/ X 

ip-^'(m,n) ip'^ '(m,n) 



(m,n) 
^2 (m,n) 



(38) 



m-1 



n-1 



m-1 



n-1 



(39) 



J=no 



We remark that we obtain the non-autonomous potential discrete modified KdV equation for ^FJI" by eliminating v™ 
from eqs.® and 

3.2 Alternate bilinearization 

There is another interesting bilinearization to the non-autonomous discrete KdV equation ([T). Let us consider the 
following bilinear equations, 



b„(am-i + a„) K^^^T^ - a„,-i(flm + b„) -d'-C^' + aja„,^i - b„) -C^/t™" 



0, 



(40) 

b„{a„,-x - flm) T^' - a„,{a,n-x - b„) t^^^i^ + a„^i(a„ - b„) t'^^i^^^^ = 0. (41) 
We obtain eq.O by introducing v™ by eq.( |35bt and eliminating i^. The Casorati determinant solution is given by 



where 



if/y{m,n) {ni,n) 



m-2 



i/Tj (in, n) 

. (s+N-X)^ N 

ij/^ (m, n) 



I {s+N-X), s 



n-X 



(42) 



i/rW(m, n) ^A.pIH + a,„p,) ]~[(1 + OiPr) ]~[(1 + bjp^) 

!=mo j=no 

m-2 n-X 
+fir(-Pry(l - OmPr) ]~[(1 - bip,-) ]~[(1 - bjp,-). 



(43) 



J=no 



Ar, fir are arbitrary constants (r - I, . . .,N) and rj" is given by eqs.(l38ll and ( [39] l. 

We note that in the autonomous case, a-" reduces to t™, the bilinear equation (l40l i yields eq.(|4]i, and eq.lHTb 
becomes trivial, respectively. Secondly, because of the symmetry with respect to m, n in eqs.([T]i and (I35bl i. the 
following bilinearization is also possible: 



a„,{bn-X + b„) C^'t™ - bn-xibn + «m) d + bn{b„-X ' aj tT^\t"Ix = 0' 

«m(/'„-i - ^) Cid' - b„(bn-x - fl,„) T;r^'c + ^«-i(^n - = o, 

whose solution is expressed as 



(44) 
(45) 



(^j''*(m, n) (^j'*^'''(m, n) 
^^(m,n) (/ij'^'^m, n) 



m-1 



(;e>j (m, n) 
(^2 {m,n) 



As+N-X). s 
0Ar '(W, «) 



«-2 



(46) 



*(^>(m,n) =i,/,;(l + Z7„p,) Y]il + aipr) ]~[(1 + bjPr) 

i=mo j=no 

m— 1 w— 2 

+flr(-pry(l - b„pr) ]~[(1 - fl/f,) ]~[(1 - V*-)" 



(47) 



]=no 



We also remark that the similar structure in the above auxiliary t functions has appeared in the study of Ri and 
Rii biorthogonal functions llT3l fT4l . Also, similar soli ton type solution has been constructed for the non-autonomous 
discrete-time Toda lattice equation lfTSl . 



We can show that tJ" and /f™ satisfy eqs.(l40]l and (l4TT l by the technique similar to that was used in refs. lfT3l[T4l[T5l . 
Namely, by using the linear relations among cp): and i^): , we first construct such difference formulas that express the 
determinants whose columns are appropriately shifted by t^' or /fjj'. Then eqs. flTb and ( l44l i are derived from Pliicker 
relations which are quadratic identities of determinants whose columns are shifted. 

From eqs.([39]l and (|43T l. we see that i^j.*' and (/r*'* satisfy 



ip^if\m + 1, n) - (fi^if\in, n) - a,„ (^''^'^(m, «), 
(p^if\m - + ^*'^'*(ot - l,n) - if^[^\m, n), 
il,\'\m, n) - fl,„ ijf^^^'Hm, «) = (!- alp]) ip\'\m - 1, «), 
<p[^\m, n+ I) - (fi'i^Xm, n) - bn ip\^*^^{m, n). 



We introduce a notation 

where 'T' denotes the column vector 



r;:'-io, 1, 



N -2, N-l 



if . (m, n) 



(48) 
(49) 
(50) 
(51) 

(52) 
(53) 



Then the following difference formulas are derived from eqs.(l48])-(l5TI) by the similar calculations to those given in 
ref.na: 





0,1,- 


■ ,N- 


3,N -2,N -2,n-i 1, 


(54) 




0,1,- 


■ ,N- 


3,A^-2,A^-2„_i 1, 


(55) 


A(m)-ifl„ t";+' = 


0,1,- 


■ ,N' 


3,A^-2,A^2„,+i |, 


(56) 


-A(m)"'(fl„_i + a,„) cr^ = 


0,1,- 


■ ,N- 


3,A?~2,„+i,A?-2,„_i 


(57) 


-(flm-i - b„-i) T^zl - 


0,1,- 


■ ,N- 


3,A?-2„_i,A?-2„_i 


(58) 




0,1,- 


■ ,N- 


3,A?-2„_i,A^2„+i 


(59) 



where 



' Ai(m)-i ^f+*\m+ 1,«) 



I A^(m)-VrV+l,n) 



A,{m)^\-aip] (r= 1,...,A?), A{m) ^\\{\ - alp].). 



(60) 



(61) 



We give the proof of the above formulas in the appendix. Applying the the difference formulas to the Pliicker relation 

0=|0,--- ,N -3,N -2,N -2„^x I x | 0, • ■ • , A? - 3, A? - 2,„_i,A?~2„,+i | 
,N -3,N -2,N -2,„^i\x\q, ~ 



0, 



,N -3,N -2,N -2„ 



X 10, 



,N-3,N-2„^uN-2„,^i 
,N-3,N-2„^,,N-2,„^i 



(62) 



we obtain the biUnear equation Equation (l44l l is derived by applying the difference formulas ( |54l i, ( |55] l, ( |58] l and 



-(am-2 - flm-l) 



-2 



-(am-2 -b„-i) cr'" I 



0, 1,-- - ,A?-3,A?-2,A?-2;„_i |, 
\0,l,--- ,N-3,N-2m^i,N^2„^i 
0, 1,-- - ,A?-3,A?-2„_i,Ar^2„,_i I 



(63) 
(64) 
(65) 



where 



V ^N^''\m - l,n) ) 



to the Pliicker relation. 



3.3 Reduction from the KP hierarchy through the potential form 

In this section, we consider the following difference equation 



1 



1 



n+l 



which is closely related to eq.([TJ as 

(1 L"! J_ - _ ,,'"+1 

\a b v'" ~ " ■ 

The autonomous version of eq.(l68]) is known as the potential discrete KdV equationfTS^. We call eq, 
autonomous potential discrete KdV equation. Casorati determinant solution to eq.(l68]) is given by 



;=mo 



J=no 



where 



Pn 



(■«)/ \ (.5+^-2)/ N (s+N), s 

ip\{m,n) ■■■ (fi^ (in,n) (m,n) 
<^2 (in,n) ^2 {m,n) 



<^2 ('W> «) 



<^)v (m,n) ip-^ '(m,n) 



and ipf\m, n)(r - 1 , . . . , A^), r™ are defined by eqs.([39]l and ( l38b . respectively. 
Equation (|68] | is derived from the following bilinear equations for p™ and t" 



m+l m _ „m _( M (^m+\ m _m+l_m^ 

\ "ah '^rt / 

„m+l /« „m / , \ 02 in+\ ni \ 

P„+l - Pn ^n+\ = —"^77 V «+l ^« ~ ^« ' 



(66) 



= 1 0, ■ 




3,N- 




|x 


o,-- 


■ ,A'- 


-3,A'- 


- 2 


i,A'-2,„-i 




-10,- 


• ,A^- 


3,N- 


2,N-2„,^i 


|x 


0,- 


• 


-3, A' 


_ 2 


i,^2,„_i 


(67) 


+ 1 0,- 


■ ,N - 


3,N- 




1 ^ 


10,- 


• 


-3,A^ 


-2„_ 


i,A?-2„,-i 


■ 



(68) 

(69) 
the non- 

(70) 



(71) 



(72) 
(73) 



through the dependent variable transformation ( iTOl i. In particular, eq.(|69]l also follows from eq.(l72]l. Therefore, we 
may regard eqs.(l72]i and (l73T l as yet another bilinearization of the non-autonomous discrete KdV equation ([T]i. 

We can show that p™ and tJ" satisfy eqs. ( l72b and ( |73] | as follows. Applying the difference formulas ( |54] |. ( fSSl l. 
and 



-(fl„_ip;r' + <-') = 10,1,- ■■ ,A^-3,A^-l,A^-2„,_i I, 
-(^_ip™, +T;f_,) = |0,l,--- ,A^-3,A^-l,A^-2„_i I, 



to the Pliicker relation 



= |0,-- - ,A^-3,A^-2,A? - 1 |x |0,-- - , - 3, - 2„-i, - 2,„-i 
-|0,-- - ,A?-3,A?-2,A^-2„_i |x |0,-- - , A^ - 3, A^ - 1, A? - 2„,_i | 
+ |0,-- - ,A?-3,A^-2,A^-2„_i |x|0,--- ,N - 3,N - l,N - 2„^i \ 



(74) 
(75) 



(76) 



we have eq. ( |72] |. Similarly, we obtain eq.(l73Tl by applying the formulas (l55T l. ( |56] |, ( |59] l. ( fTSl ) and 

A(m)-'(fl„pr' - ^r') = I 0, ■ • ■ , - 3, - 1, A^2,„+i |, (77) 

to the Pliicker relation 

= |0,-- - ,A?-3,A^-2,A?-1 |x|o,-- - , - 3, - 2„_i, A^2,„+i | 
-|0,--- ,N -3,N -2,N -2„^x I x | 0, • ■ ■ , A? - 3, A? - 1, A^2„,+i | (78) 
+ |0,-- - ,A^-3,A^-2,A^~2,„+i |x|0,-- - ,N -3,N - \,N -2„^x \ . 

We finally remark that if we introduce the continuous independent variables t\,tj„ ■ ■ ■ through ip\^\m, n) as 

m-l n— 1 

^«(m,«) ^a,pl ]~[(1 + aip,) \\{\ + bjp,.) eP'"^P'''^- 

™q J=no 

m—l n-l 



(79) 



then becomes the t function of the KdV hierarchy. In this case, p™ and m™ can be expressed as 



respectively, and m™ satisfies the potential KdV equation 

du"! 3 /<9m'«\^ 1 d^u 



P=0. (81) 



dt3 2\dhj 4 dt\ 

This is consistent with the fact that (autonomous version of) eq.(|68Tl is derived as the Backlund transformation of the 
potential KdV equation |fT6l . 

4 Concluding remarks 

In this article, we have considered the bilinearization of the non-autonomous discrete KdV equation and constructed 
Casorati determinant solution. We have presented three diff'erent bilinearizations, each of which has different origin. 
Although we have constructed only Casorati determinant solution, namely, soliton type solution, it might not be 
difficult to discuss other types of solutions, such as rational solutions or periodic solutions, based on the bilinear 
equations that have been obtained in this article. Also, we expect that other non-autonomous discrete integrable 
systems on two-dimensional lattice can be investigated in similar manner. 

As was mentioned in Section 3.2, the t functions in the second bilinearization resemble those in the theory of Rj 
and Rii biorthogonal functions, but the explicit relation is not clear yet. It might be an intriguing problem to study 
underlying structure of the second bilinearization. 

Finally, recently Takahashi and Hirota have succeeded in constructing the soliton solutions of the ultradiscrete 
KdV equation in permanent form |,17] . It might be an interesting problem to investigate the permanent type solutions 
for the non-autonomous case. 
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A Proof of difference formulas 

In the appendix, we give the proof of the difference formulas of r functions which have been used in the derivation of 
bilinear equations from the Pliicker relations for completeness. For later convenience, we first prepare the following 
two equations for i^*'' and i/'J.'' which are derived from eqs.(l48]l-(l50li: 

il/^;\m, n) + fl„,_i iAr""('M, n) = 'P^'^im +!,«), (82) 
^p['\m + 1, n) - fl„ ,^<.*+'>(m +!,«) = A,(m) ^<.*>(m, n). (83) 

Equations ([54) and ( [55) We have 

= I 0„,_i, l„,_i, ■ ■■,N- 2,„-uN- l„,_i| . (84) 
Adding the (; + l)-th column multiplied by a,„_i to the /-th column for i - 1,2, . . . ,N - 1 and using eq.(|48]), we have 

T^' =|0,1,---,A?-2,A?-1,„_,|. (85) 

Multiplying a„,_i to the A^-th column and using eq.(|48l) we obtain 

a„-i = |0,1,--- ,A?-2,fl„,_, x(A? - l)„,_i| 

^\0,l,--- ,N-2,{N-2)„-{N-2)„^i\ 
^-\0,l,--- ,N-2,N-2,„_i\, 

which is eq. (l54l) . Equation dSSl l can be shown in a similar manner by shifting n and using eq.dSTI). 
Equation ( [56) We have 

T^' = I 0„,+i, !„,+,, ■ ■ ■ , A? - 2„,+i, A? - l„,+i I . (86) 
Adding the (/ + l)-th column multiplied by to the i-th column for / = 1 , 2, . . . , - 1 and using eq. ([83] l. we have 

( Ai(m)ip^f(m,n) ^ 



=|0, 1, A^-2, A?-U+i , 



Aiv(OT) (p^^\m, n) 



Multiplying a„, to the A^-th column and using eq.(l83]l we have 

am = I 0,1, ■ • ■ ,W^,a„ x(N- 1) 



m+I 

= |0,1,-- - ,A?-2,A?-2,„+i I 
=A(ot) X I 0, 1, ■ • ■ , a? - 2, A^2,„+i I , 

which is eq.([56]l. 
Equation (|57) We have 

cr'; = |o,T,--- ,Ar^l |, (87) 
which is rewritten by using eq.([82]i from the first column to the A^-th column as 

< = I 0„,+i, • ■■ ,N- 2,„^uN^ hn I . (88) 

Now, from eqs.([50ll and ( [82] l we have 

Ar(m) ipY\m - 1 , n) = -{a„ + fl,„_i ) ip\'^^\m, n) + ipY^m +l,n). (89) 
Applying eq. ( [89] l to the A^-th column of the right hand side of eq.([88ll we obtain 
(fl,„ + fl,„_i) cr"' - - I 0,„+i, • ■ • , A^ - 2„,+i, A^ - 2„,_i I 



,N -3,n, N - 2„,+i, A? - 2„,_i I 
= -A(ot)x 0„„--- ,A^-3,„,A^-2„+i,A^-2;„_i 



which is eq. 



Equation ([58) Shifting n in eq.(l54]l we have 

- flm-iCl' = I 0„-i , l„-i , ■ ■ ■ , A? - 3„-i , N - 2„_i ,N-2 I . (90) 
EHminating (p^r^^\m - 1, n - 1) from eqs.(l48ll...-i and (BTb m-i we get 

(flm-i - b„-i) ifi';'\m - l,n - 1) = a„,-iip['\m -!,«)- b„-i(p\f\m,n - 1). (91) 
Using eq.(l9TTi to the A^-th column of the right hand side of eq.(|90l). we have 

-flm_i(fl„,_i - bn-i) t'"_i = fl,„_i I 0„_i, l„_i, ■ ■ • , - 3„-\,N - 2„-i,N - 2,„_i | . 
Adding (/ + l)-th column multiplied by fl,„_i to i-th column for / = 1 , 2, . . . , - 1 and using eq.(l48Tl. we obtain 

-(fl„,-i - b„-i) t';;:1 = i o, i, • ■ • , a? - 3, a? - 2„^i,n- 2„_i i , 

which is eq.dSSll. 

Equation (|59) Shifting n in eq.(|56]l we have 

A(m)-'fl,„ ^'^1 = I 0„-i , l„-i , • ■ • , A? - 3„-i , N - 2„_i , N^2 | 
= |0, 1,-- - ,A'-3,A^-2„_i,A^2,„+i |, 

or 

a,„ T^^l = I 0, T, • ■ • , A^, TT^n-x, N-2 I . (92) 
Eliminating ip'^^'^\m + 1, n - 1) from eqs.(l83]l„-i and (ISTTl m+\ we get 

b„-\Ar(m) ip^'\m, n - 1) + a,,, ip['\m +\,n)- (a„ + b„-i) (p'-'\m + 1, n - 1). (93) 
Applying eq.(|93Tl to the A'-th column of the right hand side of eq.(|92]i. we have 

aja,„ + bn-i) T^ll ^a,nx\0,--- ,N -3,N - 2„_i, A^ - 2,„+i | , 

which yields eq.(l59]l 

(fl„ + ^-i) T^^l = A(m) x\0,--- ,N-3,N- 2„-i , ^~2„+i | . 
We omit the proof of other difference formulas, since they can be proved in a similar manner. 
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